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Abstract 

We generalize the Toda lattice (or Toda chain) equation to the 
square lattice; i.e., we construct an integrable nonlinear equation, for 
a scalar field taking values on the square lattice and depending on 
a continuous (time) variable, characterized by an exponential law of 
interaction in both discrete directions of the square lattice. We con- 
struct the Darboux-Backliind transformations for such lattice, and the 
corresponding formulas describing their superposition. We finally use 
these Darboux-Backliind transformations to generate examples of ex- 
plicit solutions of exponential and rational type. The exponential so- 
lutions describe the evolution of one and two smooth two-dimensional 
shock waves on the square lattice. 

1 Introduction 

The Toda lattice [D EI 
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where A^/m = fm+i — fm is the difference operator and Q = Qmif) is a 
dynamical function on a one dimensional lattice (m G Z), is one of the most 
famous integrable nonlinear lattice equations. It describes the dynamics 
of a one - dimensional physical lattice whose masses are subjected to an 
interaction potential of exponential type. The infinite, finite and periodic 
Toda lattice ((T)), as well as its numerous extensions [H El El 13 El El UHl HH 
IT^ . have applications in various other physical and mathematical contexts 

o iH 113 iia im im cn • 

In this paper we introduce the following integrable generalization of the 
Toda law (Q) (i.e., the law characterized by an exponential interaction be- 
tween nearest neighbours) to a 2-dimensional lattice: 

(~1 d ( 1 dQm,-n\ _ A ((^ r< pAmQm_l,„A _|_ A ((^ p!\rQ-m,n-l\ 

''^i''^ (i%\C dt I '—^m\^in,n~^m—l,n'^ J i '—^n\^m.rv~^mp,—l'^ J ; 

(2) 

Cm+l,n+l _ ~AmA.nQm,n 
C ^ ) 

where Q = Qm,n(t) and C = Cm,n(t) are dynamical functions on the square 

lattice ((771,7?.) G Z ) and Am/m,7i ~ fm+l,n fm,ni ^nfm,n ~ /m,n+l fm,n 

are the difference operators in the m and 77 directions. 

Starting with the linear 5-point scheme (|7|) j20. , in §2 we construct a 
Lax pair for eq. (j2l); in §3 we construct the Darboux-Backlund Transfor- 
mations (DBTs) for this 2D Toda lattice, and the corresponding formulas 
describing the superposition of such DBTs; in §4 we use these transforma- 
tions to construct explicit solutions of exponential and rational type of the 
2+1 dimensional Toda lattice (0). 

We remark that, in the literature related to integrable systems, there 
exist already three 2+1 dimensional generalizations of the Toda lattice ([Q). 
The first one is obtained replacing the second derivative d? /dt^ in ^ by the 
hyperbolic operator d"^ /dxdy jH] (see also |2I1I22]), or by the elliptic operator 
d"^ /dzdz. In this equation, therefore, the scalar field Q depends on two con- 
tinuous variables x, 7/ and on one discrete variable m: Qjn{x,y). The second 
generalization [Hj can be viewed as a variant of the first, in which one of the 
two continuous variables, say x, is suitably discretized. The third generaliza- 
tion [HI [7] is obtained discretizing both x and y variables. In the generaliza- 
tion Q we propose in this paper, instead, the scalar field Q depends on the 
continuous time variable t, through the Sturm-Liouville operator in the left 
hand side of and on the two discrete indeces {m,n) G of the square 
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lattice {Q = Qm,n(t)), through the exponential law of interaction between 
nearest neighbours in both m and n directions. 

We remark that a 2+1 dimensional generalization of the Volterra system 
on the square lattice has recently appeared |23j. 



2 The 2D generalization of the Toda lattice 

The Lax pair (zero curvature representation) for the Toda lattice can be 
written in the following form 



Tm-l 

- m + 1 



^ m.-\-i ^ m i -il 



where A is the constant eigenvalue of the self-adjoint 3-point scheme 
Tmit)^ Fraif) are dynamical functions on the lattice, the eigenfunction 0m (t, A) 
solves simultaneously the Lax pair 0, and the Toda field Q is related to F, F 
in the following way 

F = r = e"^ (A) 

^ ra — (it' — ' 

A key progress towards the generalization of the Toda law to a 2- 
dimensional lattice has been recently made in [20] ; in that paper, devoted to 
the investigation of discretizations of elliptic operators on 2D lattices admit- 
ting Darboux transformations (DTs), the following results were, in particular, 
established. 

• The linear and self-adjoint 5-point scheme on the star of the square 
lattice: 

bm,n'4^m,n+l ~l~ ^m,n— 1 V'»n,n— 1 fm,n'^m,n 0, 

a natural discretization of the self-adjoint elliptic (if AB > 0) operator 

(A^^)^ + iB^y)y = (6) 

in canonical form, admits DTs. 
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(7) 



• The 5-point scheme (0) admits a distinguished gauge equivalent form: 

J- m,n+l ' i m,n ' ' ' 

obtained from (jH)) via the following gauge transformation: 

9m,n p 9m,n /„\ 

Schint — t; , [o) 

J- m,n J- m.,n 

with Fm,n = fm,n{9m,n/^m,n)^ where g and r are defined by: 

(^m,n9m+l,n ^m,n9m+l,n) ^ m,n \/ ^m,n9'm,n9iTi+l,n '\/^m,n9m,n9m,n+l • 

(9) 

The linear problem (|7j), a natural 2D generalization of the linear problem 
(jS^i.), satisfies the following basic properties: i) it possesses DTs (inherited 
from the DTs of (0)); ii) it reduces, in the natural continuous limit, to the 
2D Schrodinger equation: 

u{x, y)tp = 0. (10) 

For these two reasons the self-adjoint 5-point scheme (|7j) was identified in 
pUj as a proper " integrable" discretization of the 2D Schrodinger equation, a 
good starting point in the search for integrable discretisations of the nonlinear 
symmetries associated with the spectral problem (jlUj) and in the search for 
an integrable generalization of the Toda equation to a square lattice. 

The 2 dimensional generalization of the Lax pair Q proposed in this 
paper is indeed based on the linear problem (|7j), and reads 

r , , 'rm+l,n i p Ym—l,n ' r , , Ym,n+l i p Y^m,?i— 1 m,nYm,r, 

J- TT7,-t-1 .T7, ' J- m.n ' J- m,n+l 1 ^ ' ' ' 



~'4^m+l,n p V'm— l,n ~l~ p , , V^rriiJi+l p '^m,n~l 

It is easy to verify that this system of linear equations for the eigenfunction 
i^m,n{t) is compatible iff the coefficients r,F,C satisfy the following nonlinear 
equations: 

m+l,n \ r ,, ] '-'m— l,n I p ;~r'-^m,n+ll p I ^m,n—l\ p 

\ Am+l,n J ' \ i m,n / ' \ 1 m,ri+l / ' \ i m,n 



,71 

dT^rn.n 1 

/ \ 2 

C'm+l,n+l / '^m-\-l,n-\-l'^7n,n 



Crr 



Equations (jl2j) can be conveniently rewritten as the 2+1 dimensional gener- 
alization Q of the exponential interaction law of Toda, in terms of the scalar 
field Qm,n(t) defined by 

p ^ _ dQm,n p _ „-iQm,„ (-in) 

In the one dimensional limit in which ipm,n{t) depends trivially on n and 
the coefficients T,F,C do not depend on n: 

'^m,n(j') 0m(^)( ^) ; ^m,n Pm,n -^m; Cm; (-^^) 

it follows from (j2b) that C is an arbitrary function of t independent of m: 
Cm(t) = C(t) and equation reduces to 



1 d f 1 dQr 



C{t)dt \C{t) dt 



A^e'^"^'"-!. (15) 



Suitably rescaling t (or choosing C{t) = 1), one finally recovers the Toda 
lattice and its Lax pair Q, with A = z + z~^. 

We first remark that the 5-point schemes (jllll on the star of the square 
lattice are the simplest and natural generalization of the 3-point schemes dH}. 
Therefore we expect that our Toda 2D lattice Q be the simplest integrable 
generalization of ((T)) on the square lattice. 

We also remark that the Toda 2D-lattice system Q can be rewritten as 
a single equation, observing that equation (j2b) is identically satisfied by the 
following parametrization 

-p2 '^m,n ^ '^m-\-\,n'^m,n+\ (-\p\ 

m,n ~ ) t^m,n = 

''^r?i+l,n+l 'Tm,n'Tm+l,n+l 

in terms of the single scalar field Tm.n- Then the system (jS} reduces to the 
single equation 

'^m n+1 {j'm,r!r'm+2,n Tm+1 ,n+'fl'm-\-l ,n—l) i 

where 

W[a,l3]:=a^-a(3 (18) 

is the usual Wronskian of the two functions a and (3. It turns out j^H] that 
the scalar function Tm,n{t) is related to the r-function of the BKP hierarchy 



5 



PTj : therefore equation (|17p gives the r-function formulation of the 2+1 
dimensional Toda system Q. 

If the r-function depends only on m: Tm,nit) = Tmit), equation (fTTjl 
reduces to equation 

rl+iH[T^] - tIH[t^^,] = 0, (19) 

where 

H[Tm] ■■= frnTm - + Tm+lTm-1 - T^, (20) 

implying that 

H[rm] = fit)Tl (21) 

with fit) arbitrary function of t. By the change of variable fm{t) = 
exp{—y(t))Tm{t), with y(t) = f(t), we recover the r-function formulation 
H H[fm\ = of the Toda lattice ^. 

We finally remark that F and C constants and F = e~* is the trivial 
solution of (fT^ : correspondingly: Q = —FCt blows linearly in time. So- 
lutions which are perturbations of this solution will exhibit such a linear 
blow up in time, which can be removed introducing the change of variables 
Qm,n = Pm,n — FCt + 5. Then equation (j2)) becomes: 



-l-A (C n „A„Pm,„_i\ Cm + l,n + l _ -AmA„Pm,„ 

~r^n \_'-^m,nf-^ni,n— 1"^ ) i Cm n 

We end this section remarking that algebro-geometric solutions of the 
eigenvalue problem for a generic 5-point scheme were constructed in f2S\ . 



3 Darboux-Backliind transformations and their 
superposition 

It is a straightforward (but long) calculation to verify that the DBTs for the 
2D generalization of the Toda lattice ()12|1 . (j21) read as follows: 



Cn 



Tm — 1 ,n 
rm.n 

r 



m f) ; 



(23) 



m— l,n"m,7i— 1 
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■p/ 2 _ 

m,n 


- (r^)m-l,n-l 




c — 


(re)™_i,„(re)„ 


,n— 1 


(r9)m,n {rd)m — l 


,n — 1 ' 


F' = 

m,n 


(r^)m-l,n-l 


1 1 1 

(r6)m_l_„ (r9)m,n-l 



(24) 



+ 



\eJm-l,n ^ \ e )m,n-l 



where A_.rafm,n = fm,n-fm-i,n and A_„ = fm,n-fm,n-i- In these equations: 
dm,n{t) is a solution of the Lax pair (fTTj). for the coefficients Tm,n, Fm,n, Cm,n', 
V'm,n(^) is the transformed (via 9m,n) solution of the Lax pair (fTTj). 

for the transformed (via (|211)) coefficients L'^ ^. 

The so-called spatial part of the above DBTs was already written in 
pnj : the temporal part, describing the time dependence of the transformed 
solution ip'^^{t) of (fTTj) . and the transformation law for the coefficient 

Cm,nit), are new results of this paper. 

It is well-known j^H] that it is possible to combine DBTs of a given in- 
tegrable system, to construct superposition formulas and a permutability 
diagram of DBTs. For the 2D Toda lattice (fT^ .(j^ it is possible to prove the 
following result (see also [22] for the Bianchi permutability diagram of the 
general self-adjoint scheme on the star of the square lattice). 

Consider a solution (F^.n, -Fm.n, C*m,n) of the 2D Toda lattice (fT^ .lf^. and let 
6m,n and 6m]n bc two independent solutions of the Lax pair (fTTj) . correspond- 
ing to the coefficients Tm,n, Fm,n, Cm,n- Superimposig the two DBTs with 
respect to 6m,n and 6m,n, one obtains the new solution {Tm%, Fm^n ,Cmfn) of 
the nonlinear system ()12|) -P|) through the following formulas: 



(rm+i,ri+i)^ ~ (r 



\2 S 



m+l,n+lJ - \-m,n) ) 
p(12) _p I 1 r^-i,n //)(!) a{2) _a{l) a{2) 

1 r^,„-i fl(2) _n{l) n{2) n 

r-^ \ r?i,n— 1 m— l,n ^m—l,n m,n—l)\ 

(25) 
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where the function S^.n is obtained integrating the first order compatible 
equations 

fl ^■m-l,n^m,n-l ( oC^) f)W /d(2) N 

ujl J. m,n ' ' ' ' 

V _ V — toO^) _ fl(l) /l(2) \ (0(\\ 

V _ V — ^rn-l.n ^ gW /l(2) _ /](2) X 

This scheme is often used to construct the two-sohton solution knowing the 
one-soliton solution of the system (in this case 6m]n and 9m]n are the eigen- 
functions of the one-soliton solution corresponding to two different sets of 
parameters) . 



4 Solutions of exponential and rational type 

The existence of DBTs is considered one of the basic properties of an inte- 
grable nonlinear system. In particular, it allows one to construct iteratively 
solutions from simpler solutions, via an endless procedure. 

In this section we show some example of explicit solutions of exponential 
and rational type of the 2D Toda lattice obtained using the DBTs (j^ . 

We consider as starting solution of the system ()12j) the trivial one, corre- 
sponding to F, C constants and V = e~* (for this solution Q = —FCt) and, 
correspondingly, we look for an exponential solution of the Lax pair (jlll) : 

i^m,n{t) = e""^+^"+-*+^ (27) 

obtaing for the coefficients a,P,u! the following equations: 

cosh a + cosh /3 = -|, , . 

sinh a + sinh /3 = ^ . 

These equations can be interpreted in the following way: given the constants 
F, C, equation (j^Hk ) establishes a constraint between the "wave numbers" a 
and P; once this constraint is satisfied, equation (pSb ) gives the "dispersion 
relation" uj in terms of a and f3. 

Looking for real and non singular solutions of exponential type, in the 
following we restrict our analysis to the case F, C, a, (3,6 & M, postponing the 
study of other possible choices to a subsequent work. Then equation (|^ 
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implies that F > A and that the parameters a and P = cosh ^(f- — cosh a) 
must range in the interval 

F 

1 < cosha,cosh/5 < - - 1, F>4 (29) 

(if F = 4, then a = p = 0). 

We consider now the following solution of (jlip : 



Om,n{t) = COShe+^„(t) + pCOshe„^„(t), , . 

e^^{t) := am±f3n + iu^t + 5^, := C(sinha ± sinh/?) ^ ^ 



consisting of a suitable combination of four exponentials of the type ()27|) . 
where satisfy the constraint (j^Hk ). and p > 0, 5=*= G M. 

Applying the DBTs ()24|) to this basic solution, one obtains the following 
dressed solution of the 2D Toda lattice (fT^ .ipi: 

p/2 ^ Q-QL,At) ^ cosh e+ _i „_i (f )+p cosh e- _i „_i (f ) ^p(.^ 

cosh em,n(t)+P cosh ©m_„(t) ' f'^l 

^/ ^ [^°«he+_i Jf)+pcoshe-_i „(f)][coshe+ ,^_,(f)+pcoshe- „_i(f)] l-^^-^ 

'-^m.nl J [cosh0+,„(i)+pcoshe-,„{t)][coshe+_i_„_i(i)+pcoshe-_i_„_i(t)]^- 

This solution exhibits, depending on the value of p > 0, the following different 
features. 

i) If p = 0, we have the simplified expression: 

coshe+.„(t) ' 

^/ ^ [^°'^he+_,^jt)][coshe+,„_,(f)] l-^^J 

^m,nl J [cosh 0+,„(t)] [cosh e+_i^„_i(i)r- 

This solution describes a smooth 2D-shock wave (a kink); the shock front 
is the phase (straight) line G"*" = const, forming with the m - axis the an- 
gle —6, with 6 = tan~^^. This shock wave propagates with speed v"*" = 

— ^ sin(26')(l//3, 1/a). The values of Q'^n + ahead and behind the 
shock front are respectively —(« + /?) and a + j3\ then the shock strength is 
2(a + /?) (see Fig. 1). 

ii) If p is a finite positive number, the solution describes two smooth 2D- 
shock waves with the following features. The phase (straight) lines = 
const form with the m - axis the angles =F^; they travel with speeds = 

— ^ sin(2^^)(l//3, il/a) and, consequently, their intersection point P travels 
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with constant speed vp = — C C^™^" , ^ ) . In this situation the two shock 
fronts do not coincide, as before, with the phase hnes; they are now parallel 
to the m and n axes and intersect in P, dividing the (m, n) plane in the usual 
4 quadrants. The values of Qm,n + FCt in the first, second, third and fourth 
quadrants are, respectively, a + /3, —a + (3, — (a + and a — P (see Fig. 2). 

iii) If p is a very small positive parameter: < p << 1, the previous two 
regimes combine in the following way. In the finite (m, n) plane (or, more 
precisely, in an inner region of the order 0(ln(l/p))), the term pcoshB" is 
negligeable and the expression ()32j) is a good approximation of the solution, 
which then describes the single transversal shock wave of the regime i). In 
the outer region, that term is not negligeable anymore and the regime ii) 
becomes dominant. Both ends of the transversal shock front bifurcate into 
two semilines parallel to the m and n axes (see Fig.s 3 and 4). One could 
actually say that the regime iii) is the generic one; but, for p = 0(1), the 
inner region is not visible, since it is smaller than a single elementary square 
of the square lattice {m,n). The inner region is visible if it contains at least 
one elementary square of the lattice. If the spacing of the square lattice is 1, 
a rough extimate for this condition is that < p < mm(l/(|Q;|e), l/(|/5|e)), 
where e is the Neper constant. 

The possible existence of web-like structures in the inner region, typical of 
2+1 dimensional soliton models [30], will be explored in a subsequent work. 

Starting with the trivial solution F = 4, C constant (^ F = e^*-^*) 
of the system (fT^ . it is also possible to construct rational solutions. It is 
straightforward to verify, for instance, that 

OmAi) ■= + {Ct + a)m + {Ct + h)n + + c{a + b)t + d (33) 

is a polynomial solution of the system (jllj) . where a,b,d are arbitrary con- 
stant coefficients. Substituting it into the DBTs (j24p . one obtains a (singular) 
rational solution of the 2D-Toda lattice. We remark that this solution could 
have been derived directly from the solution for a suitable choice of its 
free parameters. 
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Figure 1: If p = and a,f3 e R {F > 4), the solution Qm,n + FCt in 
fl31|) describes a smooth 2D-shock wave propagating with velocity v+ = 
— ^ sin(2^)(l//5, 1/a). The shock front is a straight line forming the an- 
gle —9, 9 = tan~^^ with the m- axis. In this figure: a = 5, (3 = A {F = 
101), 5± = 1, C = 1, p = 



Figure 2: For a,/? G M (F > 4) and p = 0(1), the solution + FCt 
in (j31|) describes two shock waves with fronts parallel to the m and n 
axes. The intersection point P of these two fronts travels with velocity 
vp = In this figure: a = 5, (3 = A {F = 101), 5± = 

1, C = l, p = l 



13 



Figure 3: A view from the top of the solution for p = 10~^. In the central 
(finite) region, the single shock prevails; this single shock matches with the 
two orthogonal shocks, which prevail instead in the outer region. In this 
figure: a = 5, (3 = A {F = 101), 5^ = 1, C = 1, p= 10"^ 




Figure 4: A generic view of the solution for p = 10^^. In the central (finite) 
region, the single shock prevails; this single shock matches with the two 
orthogonal shocks, which prevail instead in the outer region. In this figure: 
a = 5, (3 = A{F = 101), 6^ = 1, C = 1, p = 10"^ 
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